Memory in an iV-qubit System as a Resource of the Decoherence Slow-Down with 

Application to Quantum Interferometry 
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We demonstrate that memory in an iV-qubit system subjected to decoherence, is a potential 
resource for the slow-down of the entanglement decay. As a possible application, we show that this 
effect can be used to retain the sub shot-noise sensitivity of the parameter estimation in quantum 
interferometry. 
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Introduction. — Non-classical correlations proved to be 
a key concept in many areas of modern science, such as 
quantum information 1], quantum computation, cryp- 
tography and teleportation [2] and even biology 3]. In 
ideal circumstances, non-classical states can be prepared 
and utilized at will. For instance in quantum metrol- 
ogy, correlated states can be employed to reach the Sub 
Shot-Noise (SSN) precision of parameter estimation 0- 
7]. However, a coupling with environment inevitably de- 
stroys subtle quantum correlations and drives the system 
into a classical mixture. This process, which is called 
decoherence, imposes severe limitations on what can be 
achieved in experiments [sl-fiol]. The aim of this work is 
to study the impact of decoherence on the performance 
of two-mode atomic interferometers. 

Along the line of recent interferometric experiments 
with ultra-cold atoms [lll - fl5j , we consider a collection of 
N bosonic qubits. Such system is equivalent to a single 
spin- AT/ 2 particle and therefore can be described in terms 
of the algebra of angular momentum. For simplicity, we 
consider the interferometric protocol, which consists of an 
imprint of a phase 9 between the two modes over a time 
r [lf|. In this work we take into account the coupling 
of the system to environment during the phase imprint. 
The impact of the decoherence on the preparation of the 
state has been studied in [171 ]. As our main result, we 
demonstrate - both analytically and numerically - that 
memory opens new possibilities to counteract the deco- 
herence during the evolution of the two-mode quantum 
state under the influence of the environment. 

Formulation of the problem. — A general approach for 
finding the dynamics of a system (S) interacting with the 
environment (E) relies upon solving the quantum Mas- 
ter equation [HGH, id t Qs = Tr E {[H tot , qse}}- Here we 
set H = 1 and qse is the density matrix of the system 
and the environment, while the reduced density matrix, 
Qs = TyeIqse}, is obtained by tracing out the environ- 
mental degrees of freedom. The Hamiltonian Zftot is a 
sum of three parts Hs, He and -ff; n t- The Hamiltonian 
of the system is H$ — ftJ z , where f2 is an external field 
pfjj related to the interferometric phase by 8 = Sir. Bas- 
ing on recent experimental results [H], [ijj , we assumed 
that the two-body interactions, which are essential for 
the state preparation 0, 0] , are absent during the phase 
imprint. The Hamiltonian of the environment - which 
couples with the system via H m t ~ is He- 

The procedure described above is an ab initio method, 



which in principle can always be employed to find the 
dynamics of the system. However, it is almost always im- 
practical. First of all, construction of a realistic model of 
He and H lnt is difficult because the necessary knowledge 
about the relation between numerous degrees of freedom 
is hardly accessible. Moreover, tracing out the environ- 
ment is challenging and usually gives a closed equation 
for §s only upon further radical approximations. The 
only known treatable case is the celebrated spin-boson 
model 119^, which depicts a linear interaction of a single 
qubit with a collection of independent harmonic oscilla- 
tors. 

Below we discuss a phenomenological model of the 
system-environment interaction, which circumvents the 
difficulty with determining the detailed description of the 
environment. A necessary assumption for this construc- 
tion is that the "back-action" of the system on the en- 
vironment is either absent or can be safely neglected. 
If this is the case, the only way the system can "feel" 
the environment is through an effective external field 
Q.e(£). Therefore, the interaction Hamiltonian is simply 
Hi n t — fifi(i) • J) where J is a vector of angular momen- 
tum operators. To determine the exact form of the field 
f2s(i) is as difficult as finding He- Nevertheless, statisti- 
cal properties of the environment should be much easier 
to either obtain or deduce. This observation allows for 
the final step, where we replace fi.e(f) with a fluctuating 
field described by a stochastic process u;(t) chosen to re- 
flect the aforementioned properties [2l[ . In this way, the 
explicit presence of the environmental degrees of freedom 
is mimicked by the fluctuations of the field and the total 
Hamiltonian Htot is replaced by 



H = VL.J z +u){t) J. 



(1) 



The solution of the equations of motion provided by this 
Hamiltonian gives a density matrix q, which is a stochas- 
tich process itself. It is necessary to average this result 
over the fluctuations - which is denoted by a bar - to 
obtain the density matrix of the system alone, qs = §■ 
This step is a reminiscence of tracing over the environ- 
mental degrees freedom performed in the approach of the 
Master equation. In this way, we have constructed the 
general framework of the system-environment model. To 
continue further on it is necessary to specify the statisti- 
cal properties of <*>(t). Below we argue that the station- 
ary Gaussian process is the most natural choice which 
accounts for the wide variety of physical situations. 
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Typically, we can expect that the environment consists 
of many approximately independent sources of fluctua- 
tions, which sum up to u>(t). If this is the case, than 
according to the Central Limit Theorem, u){t) tends to 
a Gaussian process with the growing number of sources. 
Therefore, the process is fully determined by its average, 
u>(t), and the correlation function, u>i(t)u)j(t') = /ey(i, t'). 
Here, subscripts i,j denote the orthogonal components of 
the vector field <*>(£). Moreover, the environment is ex- 
pected to be in a stationary state - for example in the 
thermal equilibrium. In line with the "no-back-action" 
postulate which states that the system does not dis- 
turb the environment, there is no distinguished instant 
of time. Hence, also the process itself must be station- 
ary, which implies that the average is a constant and the 
correlation function depends only on time difference. To 
make the discussion even more transparent, we choose 
u>(t) to be isotropic - which applies for an environment 
without any distinguished direction. To summarize, 



u>(t) = and Ui(t)uj(t') = n(\t - i'|)<5*. 



(2) 



where we set the average to be zero - without any loss 
of generality. In the majority of relevant situations, the 
correlation function «(|t — rapidly tends to zero when 
\t — t'\ > t c , where t c is the correlation time, which sets 
the time scale for the process [Hj]. The strength of the 
fluctuations is given by the variance of the process uIq = 

The non-zero correlation time connects the events from 
various instants of the evolution in a nontrivial way, i.e. 
not only through the initial conditions. Therefore, if the 
evolution is governed by a processes with r c ^ 0, we 
say that the system has memory. Moreover, r c can be 
considered as a measure of memory, so the larger r c , the 
more memory in the system. In this work we demonstrate 
that presence of memory is a necessary ingredient for the 
slow-down of the entanglement decay. At this point, we 
are ready to employ the above model to determine the 
evolution of the density matrix of the system in presence 
of the noise. 

Stochastic equations of motion. — The stochastic den- 
sity matrix satisfies the von Neumann equation of mo- 
tion id t g{t) = [H , g{t)], with the Hamiltonian given by 
Eq. (JXJ) . When g(t) is spanned by the appropriately cho- 
sen basis operators - which respect the symmetry of the 
Hamiltonian - the equations of motion greatly simplify. 
In our case this orthogonal basis is composed by a col- 
lection of spherical tensor operators TjfP (23| and the 
density matrix reads 



N 



m = E E 



(fk\t)Y fW 



j—0 m- 



TY T'O')tT'(i) 

_L I ± >7), ± 777, 



(3) 



where (f£\t)) = Tr f^>(t)g(0) 



The averaging over 



the fluctuations affects only (Tm\t)) on the right- hand- 
side of Eq. ([3]) , thus the problem of rinding the evolution 



of the system has been reduced to calculating these ex- 
pectation values. Although this parametrization of the 
density matrix ([3]) might seem unnecessarily complicated, 
it gives a particularly simple set of equations 



(4) 



where H^^,^) = (j, m'\H\j, m) and \ j, m) are the eigen- 
states of the angular momentum operators, which will 
be consequently used as the basis states throughout this 
work. The formal solution of Eq. ((J) is given in terms 
of the time-ordered exponential which acts on the initial 



expectation values (t$) = Tr f^ } g(Q) 



<^(*)> = E 



Texp[~iJ alt 1 H^(t') 



(5) 



For Gaussian processes, this result can be averaged an- 
alytically using the "cummulant expansion" method de- 
scribed in detail in [2JI and the outcome is 



e-m-To-O'O'+i)-"* )r+ e - 



xrmp 



(Ttf). (6) 



The decay rates are defined using a function T(t, f2) = 
Jo dti f* 1 dt 2 K(t 2 )e int2 and read T = T(t, 0), while T+ = 
Re[r(t, fi)]. The phase factor is ip = fit — T_, where 
r_ = Im[r(£, fi)]. The two rates T± result from the 
fluctuations in the x — y plane, while Tq originates only 
from the z-component of the noise. 

Equations ([3]) and ([6]) provide a general scheme for 
finding the density matrix of the system evolving under 
the stochastic Hamiltonian (JXJ) . From Eq. ^ it is evident 
that the elements of the density matrix are damped by 
the fluctuations. It will be demonstrated below, that 
the rate of damping can be substantially decreased in 
presence of memory. 

Evolution of the density matrix. — We illustrate the ef- 
fect of memory by choosing g(0) = t4')(^ , where 

I*) = 75 (If ' f ) + If ' ~f )) is the maximally entan- 
gled "NOON" state. This state decomposed into the ba- 
sis of spherical tensors in a following way 



g(0) = 



(2j) 



2 (^JV 



(AT) 



T 



(N) 
-N 



i=o Tr 



(2j)tf,(2j) 



(AT)t^(AT) 



Tr T^" n T, 



o 



(7) 



The first term is the diagonal of the density matrix 
spanned on tensors with m = 0. On the off-diagonal, 
the only non-zero elements are those in the top (spanned 
on the m = N tensors) and the bottom (m = —N) cor- 
ners. The presence of these tensors with m / in the 
density matrix of the NOON state is crucial for the SSN 
sensitivity of the J z interferometer [l6]. According to 
Eq. (J6j> the two contributions to the decay rate of the off- 
diagonal terms are the parameter-independent A^ 2 ro and 
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AT + , which depends on ft. The diagonal tensors decay 
only due to j(j + l)r + . Asymptotically only the scalar 

xf oc i remains thus the density matrix tends to a com- 
pletely mixed state j^jl which does not depend on il 
and is useless for parameter estimation. This example 
confirms that the fluctuations indeed lead to the deco- 
herence of the system. Although the decay into a com- 
pletely mixed state at large times is inevitable, we will 
demonstrate below that long memory can significantly 
slow down this process. 




FIG. 1. (color online) The temporal behavior of the damp- 
ing factor exp(— r+) for Qt c — (solid black), Q.t c = 2.5 
(dashed red) and Q.t c — 10 (dotted blue). The value of T+ is 
calculated for the Ornstein-Uhlenbeck process [27j, for which 

n(t) — ojq exp ^— Jpl-J with wot c = 1. The inset focuses on 
early times and compares the colored-noise cases with the 
white-noise limit exp (— u>oT c t) (dot-dashed green). 

In order to appreciate the role of memory as a resource 
for the decoherence slow-down, we compare the r c ^= 
case with the "white noise" limit, which is achieved by 
letting t c —> while keeping y = J °° n(s)ds constant. 
As a result, the correlation function tends to the Dirac 
delta, n{t) -> ^S(t), which gives T(t,fl) -> ^. On the 
other hand, the colored noise decay rate for short times 
1 J,J These two decay rates can be re- 



L , , 2 + 2 



uJqT c , thus A should 



is r(t,a) 

lated by noting that J^°K(s)ds 

be compared to the combination of the colored-noise pa- 
rameters uj 2 t c . Therefore, according to Eq. the white- 
noise decay is purely exponential, e~^-' +1 ^~^ , while in 
the presence of memory and for t -C t c the decay is less 

violent, namely e~^^^ +1 ^^) . As illustrated in the 
inset of Fig. [TJ which shows the temporal behavior of 
the damping factor at short times, the memory unlocks 
a new time-scale for experiments at which the decoher- 
ence build-up is slow. In contrast, the memory-less white 
noise case is characterized by the omnipresent exponen- 
tial decay of coherence. 

Another way to control the decoherence results from an 
interplay between the deterministic field and the fluc- 
tuating part u){t). To picture this effect, we note that 



every realization of u(t) yields a realization of stochastic 
g(t) as a solution of the von Neumann equation of motion. 
Each realization can be viewed in an approximate "stro- 
boscopic" picture, where the field u>(t) does not change 
much within correlation time r c , and then jumps to an- 
other value. The aforementioned interplay becomes clear 
when the evolution is viewed from the reference frame 
rotating around Q. Transition to this frame does not 
change the z component of the fluctuations which is man- 
ifested by the f2- independent decay rate r . On the 
other hand, the perpendicular part rotates around the 
z-axis with frequency Q. If the period of the rotation is 
small in comparison to the correlation time (£It c 3> 1), 
it performs many revolutions between jumps. In such a 
case, stochastic jumps are not noticable anymore, differ- 
ent realizations become indistinguishable thus one can- 
not speak of fluctuations anymore. This gyroscopic effect 
is exhibited by the decrease of r+, as shown in Fig. [TJ 
which compares the damping factors e~ r + for three dif- 
ferent values of Ot c . Note that in the white-noise limit, 
the jumps are instantaneous which leaves no time for the 
gyroscopic effect to kick in and consequently the decay 
rate A does not depend on fi. 

In the next step we examine how the build-up and 
gyroscopic effects can be utilized to preserve the SSN 
sensitivity in the interferometric parameter estimation. 

Memory in noisy quantum interferometry — The preci- 
sion of the estimation of the parameter from a series 
of n experiments is bounded by the Cramer-Rao Lower 
Bound (CRLB) [H, 



A 2 n> -— 



nF, 



Q 



(8) 



where Fq is called the Quantum Fisher Information 
(QFI), which depends on the state q(t) used for esti- 
mation, namely 



(9) 



Here, \i) denotes the eigen-state of q(t) with a corre- 
sponding eigen- value pi . When Fq > Nt 2 , the two-mode 
state q{t) is usefully entangled from the interferometric 
point of view [!, 0] and via the CRLB ([8]) allows for the 
SSN sensitivity. In the absence of noise, it is the NOON 
state used in the previous section that provides the max- 
imal estimation precision, as Fq — N 2 t 2 reaches the 
Heisenberg Limit (HL) 0, 0j- As argued above, when 
this state is exposed to noise for a long time, it reaches a 
completely mixed state, which gives Fq = and thus is 
utterly useless for interferometry. However, in presence 
of memory effects, so that NT + can be considered small, 
the N 2 scaling of the QFI which reads 



N 2 



1 - AT 



| 2 e -2jv 2 r -2ivr% 



+ 



NT+ [d n k^ 



(10a) 
(10b) 
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decays slower than in the white-noise limit 

F ^wn) =iv 2 T 2 e -2JV(JV + l)^ > (n) 

Moreover, additional information about f2 comes from 
the dependence of the decay rates on the parameter. We 
stress that this improvement is purely a memory effect, 
absent in Eq. ([TT|) . In both cases, the strongest source of 
decohcrence - which cannot be diminished by the gyro- 
scopic effect - comes from the z-component of the noise 
ui z (t) which sets Tq. Note that in the Hamiltonian ((T|), 
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FIG. 3. (color online) Illustration of the structure of the ini- 
tial density matrices of the states generated with the Bose- 
Hubbard Hamiltonian. The contour plot shows the distri- 
bution £\ |Tr (e(0)f, ( „ j) ) | 2 /Tr (Ti 3)t ri j) ). For y ~ 0, the 

density matrix consist of close-to-diagonal terms, so only the 
tensors with m ~ contribute. When 7 < —1, there is a 
growing contribution from large m's. 



FIG. 2. (color online) The QFI in units of r 2 for different 
input states g(0) of N — 100 particles parametrized with 7 
of the Bose-Hubbard Hamiltonian (see text for details) and 
with u) z (t) = 0. The Figure compares the QFI at r = 0.5r c in 
absence of noise (solid black), in presence of memory (dashed 
red) and in the white- noise limit (dotted blue). The grey 
horizontal dot-dashed line denotes the shot-noise limit Fq — 
100. The inset shows the QFI for the isotropic case in presence 
of memory (dashed red) and in the white-noise limit (dotted 
blue). As in Fig.[T]the colored noise is the Ornstein-Uhlenbeck 
process with ljot c = 1 and Qt c = 10. 

to z (t) can be interpreted as fluctuations of the estimated 
parameter itself. In this spirit we relax the isotropic as- 
sumption {2J and let the strength and correlation time 
of cj z (t) be different from the perpendicular components. 
It is reasonable to assume that the experiment is set up 
in such a way, that the estimated parameter is well de- 
fined. That is, the fluctuations along the z-axis should 
be slow and weak. In particular, by setting uj z = 0, we 
get To = so the strongly damping exponential factor in 
line (|10a|) disappears. Also, in the white-noise limit the 
decay rate is less violent and equal to 2N/T rather then 
2N(N + 1) /T in Eq. (fTT]). 

To picture the impact of noise, in Fig. [2] we plot the 
QFI for a wide family of usefully entangled input states 
£?(0). These states are the ground states of the two- well 
Bose-Hubbard Hamiltonian Hbh = —EjJ x + UJ 2 with 
attractive two-body interactions jiH [26j j for different val- 
ues of the ratio 7 = U/NEj. And so, for 7 = we obtain 
a coherent state, while for 7 — > —00 it tends to a NOON 



state. The structure of this family of states in terms 
of the spherical tensors is presented in Fig. [3J Figure [5] 
shows that for (J z (t) — 0, memory effects prevent the sig- 
nificant information loss and keep the particles entangled, 
as Fq > Nt 2 . To contrary, in the white-noise case, the 
SSN sensitivity is lost even at short times. Moreover, the 
inset underlines the destructive role of the fluctuations of 
the parameter f2. 

Conclusions and acknowledgements — We have demon- 
strated how the memory in an A^-qubit system, which is 
a subject to external noise, can be employed to suppress 
the entanglement decay. We have argued that this ef- 
fect has an application to quantum interferometry. In 
presence of memory, the deterministic evolution - which 
imprints the parameter- dependence on the state - cou- 
ples to the non-unitary evolution generated by the fluc- 
tuations. As a result, the SSN sensitivity is retained 
for longer times - not only due to the decoherence slow- 
down but also thanks to additional information about f2 
present in the decay rates. 
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